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Abstract 

We present a calculation of the meson masses, decay constants and quark-antiquark 
vacuum expectation value for the three generic QCD-like chiral symmetry breaking 
patterns SU{Nf) x SU{Nf) SU{Nf)v, SU{Nf) SO{Nf) and SU{2Nf) 
Sp{2Nf) in the effective field theory for these cases. We extend the previous two-loop 
work to include effects of partial quenching and finite volume. 

The calculation has been performed using the quark flow technique. We reproduce 
the known infinite volume results in the unquenched case. The analytical results can 
be found in the supplementary material. 

Some examples of numerical results are given. The numerical programs for all 
cases are included in version 0.54 of the CHIRON package. 

The purpose of this work is the use in lattice extrapolations to zero mass for 
QCD-like and strongly interacting Higgs sector lattice calculations. 


1 Introduction 


Effective field theory is used extensively in the study of strongly interacting gauge theories. 
A recent review covering a number of different applications in addition to other methods 
is [I]. Besides general interest in understanding strongly interacting gauge theories, they 
might still be useful as an alternative for the Standard Model Higgs sector as well as for 
dark matter. These applications have been reviewed recently at the 2015 |21 E] and 2013 
[1] lattice conferences. A number of recent lattice studies is [5] . Reviews of technicolor and 
strongly interacting Higgs sectors are [SlCli- 

Lattice studies are always performed at a nonzero fermion mass. In order to obtain 
results in the massless limit extrapolations are needed. A main tool for this in the context 
of lattice QCD is Chiral Perturbation Theory (ChPT) P [TII| ITT]. 

In the case of equal mass fermions three main symmetry breaking patterns are possible 
[la [m [n]. For Np Dirac fermions in a complex representation the global symmetry 
group is SU{Nf)l x SU{Nf)r and it breaks spontaneously to the diagonal subgroup 
SU{Nf)v- For Dirac fermions in a real representation the global symmetry group is 
SU{2Nf) and it breaks spontaneously to SO{2Nf). An alternative possibility is that we 
have Np Majorana fermions in a real representation with a global symmetry group SU{Nf) 
spontaneously broken to SO{Nf)- We show in this work that the EFT for the quantities 
we consider is really the same as for Dirac fermions. The hnal case is Np Dirac fermions 
in a pseudo-real representation. The global symmetry group is again SU{2Nf) but in this 
case it is expected to be broken spontaneously to Sp{2Np). 

The effective held theory (EFT) for these cases is discussed at tree level or lowest order 
(LO) in [15]. At next-to-leading order (NLO) the hrst case is simply ChPT for Np light 
quarks with a symmetry breaking pattern of SU^Np) x SU{Nf) —)■ SU^Np), a direct 
extension of the QCD case and was already done in The pseudo-real case was done 
at NLO by [IB]. The SU{2Nf) —)■ SO{2Nf) case was done in [T7]. The extension for all 
three cases to next-to-next-to-leading order (NNLO) was done in earlier work by one of 
the authors [18]. More references to earlier work can be found there and in [191 EQ]- 

This paper is an extension to the work of [H]. We add a short discussion showing 
that the calculations and the Lagrangian for the real case also covers the case of Majorana 
fermions in a real representation. The main part of the work concerns the extension of the 
calculations at NNLO order of the masses, decay constants and vacuum expectation values 
to include effects of partial quenching and hnite volume. 

Partial quenching was introduced in ChPT by 121 . A thorough discussion of the as¬ 
sumptions involved is in [22]. It allows to study a number of variations of input parameters 
at reduced cost, as discussed in e.g. [23]. We do not use the supersymmetric method in¬ 
troduced in [2T] and extended (at NLO) to the cases discussed here in [T7]. We only use 
the quark-flow technique introduced in m- Two-loop results in inhnite volume partially 
quenched ChPT (PQChPT) for the masses and decay constants are in [25l 1261127] . The 
dehnitions of the inhnite volume integrals we use can be found there. 

Finite volume effects in ChPT were introduced in ChPT in [25112B| I5U] . Early two- 
loop work is [SHEg. The vacuum expectation value was discussed in more detail in [33] . 
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After the proper evaluation of the hnite volume two-loop sunsetintegrals using two different 
methods [M] the masses and decay constants were treated in both the unquenched [3S] 
and partially quenched [36] case. In particular the integral notation at hnite volume we 
use is dehned in [36]. 

In Sect. [2] we recapitulate briehy the discussion from [18] at the quark level and add 
the case with Majorana fermions. Sect. [3] similarly recapitulates [13] at the effective held 
theory level and adds the Majorana fermion case. The cases with Dirac fermions and 
Majorana fermions are essentially identical from the EFT point of view for the quantities 
we consider. The underlying reason is an U{2Nf) transformation that relates the two cases 
as discussed in Sect. ID Partial quenching and the quark.how techniques we have used for 
the diherent cases is discussed to some extent in Sect. [5l For a discussion on hnite volume 
and the notation used there we refer to [33] • Our analytical results are described in Sect. E] 
in particular we clarify the dehnitions of the decay constant and vacuum expectation value 
used in terms of quark helds. The numerical examples and checks are presented in Sect. [71 
The analytical formulas are included in the supplementary hie 1371 and the numerical 
programs are available via CHIRON, [3HI [39]. The last section briehy recapitulates the 
main points of our work. 


2 Quark level 

2.1 The three Dirac fermion cases 

The discussion here is kept very short, longer versions can be found in [15] and [18]. This 
subsection is mainly included to show normalization conventions. 

QCD or complex representation In the Np equal mass Dirac fermions in a complex 
representation, we put the Np fermions together in an Np column matrix q. The global 
symmetry transformation by Qp X SU{Np)p x SU{Np)ji is given by 

Qr ^ ^ 9Llii9L + i9LdiMgL , ^ gplfj^gh^ + igpd/j^g]^, Ai^gnAigi- 

( 1 ) 

The matrix At = rugl + s + ip brings the quark mass term rriql and the external scalar 
s and pseudo-scalar densities in the Lagrangian via —qpM.qL + h.c.. The external helds 
are in the Lagrangian via qp'y^l^qL + QRl^'^^fMQR- Taking derivatives w.r.t. the external 
helds allows to calculate relevant Green functions mm- In particular, deriving w.r.t. Sn 
allows us to obtain {qpiqRi + qRiqLi) and derivatives w.r.t. a^i 2 with r^i 2 = —lfj .12 = 0,^12 
allows access to matrix-elements of ^ 27 ^ 7 ^*?! The symmetry is spontaneously broken by a 
vacuum expectation value 

{liflR-^ = ■ ( 2 ) 

This leaves a global symmetry SU{Np)v with gp = gp unbroken. 
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Adjoint or real representation When the fermions are in a real representation, we can 
introduce besides the Np right handed fermions qm a second set of right handed fermions 
in the same gauge group representation, qm = CqJ^p These can be put together in a 2Nf 
column vector q, q^ = {qm ... qpNp qm ■ ■ ■ qRNp)- The global symmetry transformation 
with g G SU{2Nf) is now 

Q -t gq, 9 V^,g\ M gMg'^. (3) 


We dehne the external densities and currents as in the QCD case with and Ai. We 
dehne 2Nf x 2Nf matrices 


M 


0 \ 

0 J’ 




(4) 


Note that the global symmetry can change quark-antiquark currents to diquark currents. 
The fermions condense forming a vacuum expectation value 


1 

2 


(qfc’qj) 


wTsp -Js — 



This leaves a global symmetry SO{2Nf) with gJsg^ = h 


(5) 


Ac = 2 or pseudo-real representation When the fermions are in a pseudo-real rep¬ 
resentation, we can introduce besides the Np right handed fermions q^a again a second 
set of right handed fermions in the same gauge group representation, qjna = ^abCqJ^ip a, b 
are gauge indices and the extra Levi-Civita tensor eab is needed to have qpia transform 
under the gauge group as qipa- The explicit formula is for the case of the fundamnetal 
representation with Nc = 2. qpi and qm can be put together in a 2Np column vector q, 
q^ = {qm ■ ■ ■ qRNp qRi ■ ■ ■ qRNp)- The global symmetry transformation with g G SU{2Np) 
is now 

q gq, % g%g\ M gMg^- ( 6 ) 


We dehne the external densities and currents as in the QCD case with r^,/^ and A4. We 
then dehne 


M = 




(7) 


Note that the global symmetry can again change quark-antiquark currents to diquark 
currents. The fermions condense forming a vacuum expectation value 


2^{qja^abC'qjb'} J A 


0 

I 



This leaves a global symmetry Sp{2Np) with gJAg"’" = h 


( 8 ) 
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2.2 Majorana fermions in a real representation 

In the earlier work [18] at infinite volume Dirac fermions and Dirac masses were assumed. 
It was then also asumed that the vacuum condensate was aligned with the Dirac fermion 
masses. There is in fact another possibility. Majorana fermions with a Majorana mass in a 
real representation of the gauge group. In this case the global symmetry is SU{Np). It is 
expected to be spontaneously broken down to SO{Np) which is aligned with the Majorana 
masses. 

A Majorana spinor is a Dirac spinor that satisfies 

i, = ct or ) , (9) 

The last equality are in the chiral representation for the Dirac matrices. The Lagrangian 
for a single free Majorana fermion is 

^ . ( 10 ) 

= I, = —ah If we want to gauge this for m 7 ^ 0 the mass term requires the fermions 
to be in a real representation of the gauge group. 

For Np Majorana fermions V’mi in fhe adjoint representation with external fields 1^ 
and A4 the Lagrangian, put in a big column vector q'^ = is 

C = - ^tr^ + gV^ATg*) . (11) 

This Lagrangian has a global SU{Np) symmetry with g G SU{Np) with 

Q ^ 9Q,, 9y^l9^ + igdi^g ^, M gMg^ ■ (12) 

The maximal symmetry argument says that in this case the fermions will condense to 
the flavour neutral vacuum (trc (g^Cg)). This is conserved by the part of the global group 
that satisfies gg^ = I or the conserved part of the global symmetry group is SO{Np). 

Note that the form of the vacuum and the form of the mass term are the only differences 
as far as the global symmetry group and its breaking are concerned compared to the case 
with Np/2 Dirac fermions in a real representation. 

3 Effective field theory 

3.1 The general LO and NLO Lagrangian 

The ChPT Lagrangian for Np flavours at LO and NLO has been derived in m- The 
Lagrangian for the other cases has the same form as has been shown in [151 HH] and other 
papers. The precise derivation can be found in [181 and the Maiorana fermion case below 
in Sect. ESI 


4 



In terms of the quantities f±^v^ x± defined below for each case the lowest order 
Lagangian is 


C2 = ^{u^u^+ X+) ■ (13) 

Here we use the notation (H) = tri7’(H), denoting the trace over flavours. The NLO 
Lagrangian derived by HU reads 


£4 = Lq (u^u’^u^u^) + Li {u>"u^) {u^Uy) + L2 {u^u'') {u^Uy) + Lg {u^u^u''uy) 

+La{u^u^){x+) + L^{u^u^x+) + Lq(x+)‘^ + Lt(x-)‘^ + -^L^{xX + X-) 

—'iLg{f+^_iuU^u'^) + ^Tio(/+ — fl) + — /^) + -H2{xX ~ T-) • (14) 

The NNLO Lagrangian has been classihed for the A^iT’-flavour case in M- The Lagrangian 
at NNLO for the other cases is not known, the direct equivalent of the results in M is 
dehnitely a complete Lagrangian but might not be minimal. For this reason we do not 
quote the dependence on the NNLO Lagrangian in the real and pseudo-real cases. 

The divergences at NLO were derived for the QCD case in [lU , for the others in [TBlfTS] . 
At NNLO only the QCD case is known [TT] . 


3.2 The three Dirac fermion cases 

A more extensive discussion can be found in HSIIIB]. Here we simply quote the results. 

When we have a global symmetry group G with generators T“ which is spontaneously 
broken down to a subgroup H with generators which form a subset of the T“, the 
Goldstone bosons can be described by the coset G/H. This coset can be parametrized [12] 
via the broken generators Below we explain what is used for the different cases. We 
always work with generators normalized to 1, i.e. = 6°''^. 

The quantities used from the quark level are given in Sect. [2j 


QCD or complex representation The Goldstone boson manifold is in this case SU{Np) x 
SU{Nf)/SU{Nf) which itself has the structure of an SU{Nf) Note that the axial gener¬ 
ators do not generate a subgroup of SU{Nf) x SU{Nf) even if G/H has the structure of 
a group in this case. 

We choose as the broken generators the generators of SU{Nf) ~ G/H. The 
quantities needed to construct the Lagrangian and their symmetry transformations are 


u 

Uf, 

X 

X± 


“ATS? 

i - iri,)u - u{df, - Q)m^) -)■ hUf,h ^, 

2BqM. —)■ QRXdL 

± ux^u —)■ hx±h^ , 


^ -)■ gnuh"^ = hugl 
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rf^u 

f±fip 

The first line defines h 


dfjlu HimIim T ilulfi ^ QL^fiuQL 

df,r^ - df,r^ - ir^Vf, + 

± u'r^yu -)■ hf±^^h'' . 


( 15 ) 


Adjoint or real representation The Goldstone boson manifold is in this case SU{2Np)/SO{2Np). 
The unbroken generators satisfy Q°^Js = which follows from gJsg^ = Js- The 

broken generators satisfy JsX°' = X°‘'^ Js- 

The quantities needed to construct the Lagrangians are [IB] 

u = 

X = 

X± = 

V^U = 

The hrst line dehnes h by requiring that guh^ is of the form exp(z7r“X“/(\/2F)). Note that 
the derivation used Jpu = Js- 


i {u\d^ - iVf,)u - u{d^ + , 

2BoM 

u^xJsu^ ± uJsx'^u 

- t - VuV^?j 

JsuVf.^uUs ± (16) 


^ —)■ guh"^ 


Ac = 2 or pseudo-real representation The Goldstone boson manifold is SU{2Np) / Sp{2Np)- 
The unbroken generators satisfy Q^Ja = —JaQ°^^ which follows from gJAg^ = Ja- The 
broken generators satisfy JaX°^ = X°''^Ja- 
The quantities needed are [T8] 


u 

X 

X± 

V^.U 

f±liu 


exp ( ^ 7r^X° 

\V2F 

i (u^idf, - iVf,)u - u{di, + iJAV^J aW) ^ 

2BoM 

u'xJ'a'^^ T uJa')^u 

d,Vu - d^V, - t (l> K - KV,) 

JaFVJ i UV^l/U^ 


^ —)■ guh^ 


(17) 


The hrst line dehnes h by requiring that guh^ is of the form exp(i7r“A“/ (\/2F))- Note that 
the derivation used Jau = Ja- 
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3.3 Majorana fermions in a real representation 

The vacuum in this case is characterized by the condensate 

\{QlCqj) = ^{qq)Sij . (18) 

Under the symmetry group g G SU{Nf) this moves around as 

Sij { 9 ^ 9 ).. ■ (19) 

The unbroken part of the group is given by the generators and the broken part by the 
generators which satisfy 

Qa ^ _QaT ^ ^ j^aT _ ^ 20 ) 

Just as in the cases discussed in [18] we can construct a rotated vacuum in general by using 
the broken part of the symmetry group on the vacuum. This leads to a matrix 


U = uu^ —)■ gUg^ with u = exp 

The matrix u transforms as in the general CCWZ case as 

u —)■ guh^. 


.\/2F 


TT A 


( 21 ) 


( 22 ) 


Some earlier work used the matrix U to describe the Lagrangian [13]. Here we will, as 
in [TS] use the CCWZ scheme to obtain a notation that is formally identical to the QCD 
case. We add Np x Np matrices of external fields and A4. We need to obtain the u^, or 
broken generator, parts of — iV^) u. Eq. fl^ have as a consequence that u satisfies 


u = u 


(23) 


This leads using the same method as in [18] to 

= i - iV^)u - u{di, + . (24) 

With this we can construct Lagrangians. The equivalent quantities to the field strengths 
are 

f±^iv = ± uV^uv) (25) 

with = dfjYiy — dyV^ — i {v^Vy — VyVi^j and for the mass matrix 

X± = ± u^x'u (26) 

with y; = 2BqM.. The Lagrangians at LO and NLO have exactly the same form as given 
in flT^ and (ITT]) with x± ci-nd as defined in (12^ and (126]) . 
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4 Relation Dirac and Majorana for the adjoint case 


As discussed below, we have calculated the adjoint case using two methods. They were 
appropriate for the Dirac and the Majorana case respectively. After doing the trivial 
2Np Np change the results agreed exactly. If we compare the two cases, we see that 
the main difference is really the choice of vacuum. 

The Dirac and Majorana cases lead to a choice of vacuum 

{qIC qj )DO^Jsij, {qf Cqj )d oc 1^ . (27) 


Is it possible to relate the two cases in a simple way? Under a global symmetry transfor¬ 
mation the hrst one transforms as Js —)■ gJsd"’"■ If we could hnd a global transformation 
gp that lead to gpJsgp = I the two cases would be obviously the same. 

It is not possible in general with a SU{2Np) rotation to accomplish this since det Js = 
±1 (—1 for the 2Np = 2) while det I = 1. However it is possible with a U{2N) transfor¬ 
mation. An explicit choice for gp, with a free phase a is 


_ Tie*"I ±ie-*“I \ 

9k - e'"I e-“I J 


(28) 


It can be checked that this transforms a Dirac mass term for Np Dirac fermions into a 
Majorana mass term for 2Np Majorana fermions. 

Inspections of the effective Lagrangians needed lead to the immediate conclusion that 
the mass independent terms really are U{2Np) invariant, and the mass dependent terms 
for the two cases are turned into each other. 

gp can also be used to relate the two different embeddings of SO{2Np) in SU{2Np) to 
each other. For the Dirac case the SO{2N) generators satished Q°'"’"Js = —JsQ°' while for 
the Majorana case they satished = —Q°‘. The two sets of generators are related by 

< 3 “ = gH 0 “gli. .v = gHA'“9P (29) 


5 Partially quenching and the quark flow technique 

A thorough discussion of PQChPT and in particular the derivation of the propagator used 
there is [13]. That discussion uses the supersymmetric method. Alternative methods of 
calculation are the replica trick jH] and the quark how method [21] . The earliest partially 
quenched work for QCDlike theories used the supersymmetric method DU- The replica 
trick has been used in [15]. We use the quark-how method. 

For this method we look at the matrix 

$ = (30) 


for each of the cases. 

For the QCD case, $ is a traceless Hermitian matrix. We actually keep <I> in the havour 
basis with elements (pij and i, j are havour indices. The tracelessness condition is enforced 


by the propagator. The indices are kept explicitly and the propagator connecting a held 
(pij to (pki is |13] 

Gijkiik) = Gij{k)6ii6jk - 6ij6kiGj^{k)/rises, ■ (31) 

The number of sea quarks risea is what we call Np. with G^j = i/{p^ — Xij)- The neutral 
part of the propagator, Gjj^, can contain double poles. In particular for the mass cases we 
consider: 


Glp = i{Xi-X4)/{p^-Xif+ i/{p^-Xi), 

Gls = V(P^-Xi), 

Gls' = V(P^-X4)- (32) 

n, s denote valence or sea quarks. The extra parts come from integrating out the <I>o |13] and 
enforce the condition that $ must be traceless. When constructing the Feynman diagrams, 
we keep all Favour indices free. Those that connect to external states get replaced by the 
value of the external valence Favour index and the remaining ones are summed over the 
sea quark Favours. In the present calculation, with all sea quarks the same mass, that 
corresponds to a factor of Np for each free Favour index. 

For the Majorana, SU{Np) —)■ SO{Np), case we have that $ = with $ Hermitian, 
traceless and symmetric. Hermitian and traceless follow from SU{Np) and symmetric from 
(120]). Going to the Favour basis for the diagonal elements of $ there is no change w.r.t. the 
QCD case, but the Favour charged or oF-diagonal elements must be correctly symmetrized. 
This has to be done both for the propagator and the connection to the external states, 
keeping track of the needed normalization. Afterwards we set the Favour indices connected 
to external states to their valence values and sum over the Favours for the free indices. 

For the Dirac adjoint case, SU{2Np) —)■ SO(2Np), case we have that <h = 7r“X“ with 
<F Hermitian, traceless and satisfying X°‘Js = and the matrix $ is 2Np x 2Np. 

Rewriting d) with Np x Np matrices leads to the form 

4>=(|j I?), with {/1>=0, 4c = 4j. (33) 

is Hermitian. The elements in correspond to quark-antiquark states, those in <Fc to 
diquark states. can be treated exactly as in the QCD case, both the diagonal and Favour 
charged or oFdiagonal elements, since (‘I’a) = 0 replaces (d*) = 0 in the QCD case. <I)c 
can be treated as oFdiagonal or Favour charged propagators but the needed symmetrizing 
should be taken care of both for external states and propagators. The normalization of 
all states must be done correctly as well. After constructing Feynman diagrams with both 
<h/i and <I>c degrees of freedom taken into account, we sum free index lines over the Np 
degrees of freedom, not 2Np. The results always agree with the calculations done with the 
previous, Majorana, method. 

For the last case, SU{2Np) —)■ Sp{2Np), pseudo-real, we have that d* = with 

$ Hermitian, traceless and satisfying and the matrix d) is 2Np x 2Np. 
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Rewriting $ with Np x Np matrices leads to the form 



with (A) = 0 , <hc:> = —< 1 *^ . 


(34) 


is Hermitian. The elements in correspond to quark-antiquark states, those in $<7 
to diquark states. 3)^ can be treated exactly as in the QCD case, both the diagonal and 
flavour charged or offdiagonal elements, since (^a) = 0 replaces (<h) = 0 in the QCD 
case. <hc can be treated as offdiagonal or flavour charged propagators but the needed 
antisymmetrizing should be taken care of. The normalization of all states must be done 
correctly as well. After constructing Feynman diagrams with both and $ 0 - degrees 
of freedom taken into account, we sum free index lines over the Np degrees of freedom, 
not 2Np. In this case and the previous we can also compare calculations with $a or 
external states providing a check on our results. 

6 Analytical results 

We have calculated the masses, decay constants and vacuum expectation values to NNLO 
for the QCD-like theories with the symmetry breaking patterns discussed above. A number 
of checks have been performed on the analytical formulas. The inhnite volume unquenched 
results were obtained earlier in [TS] and we have reproduced those. The partially quenched 
and hnite volume results in the QCD case are hnite. The partially quenched expressions 
reduce to the unquenched results whenever we set the sea mass equal to the valence mass. 
In addition we reproduce the known results at NLO for the condensate UTl also for the 
partially quenched case. The hnite volume expressions have been checked against the 
known NLO results and numerically with the earlier known NNLO results, as discussed in 
Sect. [71 

For the real and pseudo-real case we have the additional check that calculating the 
mass or decay constant of a quark-anti-quark or a diquark meson gives the same results. 
This corresponds to using a held from the A or the C sector in the matrices fl33ll34p . For 
the real case we have the additional check that the results using the Dirac case and the 
Majorana case coincide. 

The hnite volume case is always done for three spatial dimensions of size L and an 
inhnite temporal volume. In addition we work in the center of mass system, the momenta 
are such that the external states have zero spatial momentum. 

The masses are the physical masses as dehned as the pole of the full propagator. We 
consider here the case where all valence quarks have the same quark mass mi = rh and the 
sea quarks all have the same mass = mg- For the unquenced case obviously = mi. 
The labeling is similar to those used in three havour PQChPT [25l [26l ETJ [36]. In the 
formulas we use instead the quantities 


Xi = 2Ro"ii, X4 = 2 RoX4, Xi 4 = - (Xi + Xa) ■ 


(35) 
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These quantities are referred to in [27] as mil, m44 and ml4 respectively. 

The formulas are given for the cases SU{Np) x SU{Np) —)■ SU{Np), SU{Np) 
SO{Np) and SU{2Np) —)■ Sp{2Np). Note the difference in convention for the second case 
compared to [Hj. The three cases are referred to in the formulas with SUN, SON and SPN 
for the unquenced case and PQSUN,PQSON and PQSPN for the partially quenched case. 
In the latter case Np referes to the number of sea quarks. 

For the mass we consider a meson made of a different quark and anti-quark or a diquark 
state with two different quarks. These are always valence quarks. The physical mass at 
hnite volume is given by 

rUphys = Xi + . (36) 

The superscript (n) labels the order correction and A'^ indicates the hnite volume 
corrections. In all cases the lowest order mass squared is given by xi. A further break up 
is done for the LEG dependent parts via the (NLO) and (NNLO) and the remainder 
via 


^(4)2 ^ 

^(6)2 ^ ^i^(6)2^^L(6)2^^iJ(6)2 

All quantities are given explicitly in 133. 

The decay constant Aphys for the same mesons as above is expanded w.r.t 
order as 

^Thys = Flo (l + F^^^ + , 

with a similar split 

_p(4) _ pL{4) _|_ pR{4) 

_p(6) _ pK{6) pL{G) _|_ pR{6) 

^VpiG) ^ ^Vpm ^ ^VpR{6) (^32^ 

All quantities are given explicitly in 133- 

The vacuum expectation value is expanded in exactly the same way 

l^phys = VlO (l + -f , (40) 

with a similar split 

^(4) ^ ympyRi^) 

y(6) ^ ^K{6) _j_ ^L{4) ^R{4) 

( 41 ) 

All quantities are given explicitly in [37] . 


. the lowest 
(38) 
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The quantities with K for the SON and SPN case have been set to zero. They are 
polynomials up to the needed degree in xi and X 4 , with an overall factor of Xi for fhe 
mass. 

The decay constant and the vacuum expectation value were dehned implicitly in [18] 
using a generator X°‘ in the axial current normalized to one and an element in A4 normal¬ 
ized to one. The consequence was that in [18] Tlo = F and vlq = —BqF^ for all cases. 
This is not exactly what was done in earlier work leading to differences in factors of 2 and 
y/2. Below we explicitly specify all dehnitions in terms of the quark helds. 

QCD or complex representation If we label the hrst Dirac (valence) quark by 1 and 
the second by 2 the decay constant and vacuum expectation value are dehned as 

(0|gi7/^75<?2|Tf(p)) = iV2FphysPf, 

{QiQi) = {QliQri + QriQli) =Vp^ys (42) 

M denotes a meson of that quark content with momentum p. 

The resulting lowest orders are 

Flo = F vlo = —BqF^ . (43) 

Adjoint or real representation Here we have to be careful how we dehne the physical 
decay constant. We can choose to do using generators normalized to one using Dirac 
Fermions or generators normalized to one using the g* elements. 

With a Dirac fermion dehnition, the hrst Dirac (valence) quark labeled by 1 and the 
second by 2 , the dehnitions are 

(0|gi7/.75g2|M(p)) = i\/2FphysP;, 

{QiQi) = {QliQri + QriQli) =Vpi,y^ (44) 

M denotes a meson of that quark content with momentum p. The resulting lowest orders 

are 

Flo = V2F vlo = -2BoF^ . (45) 

If we instead choose to use the Majorana case, the natural dehnition of the decay 
constant and vacuum expectation value with the hrst (valence) Majorana fermion labeled 
as 1 and the second as 2 via 

■^^{0\Q*lF^q2 + q2^^,ql\M{p)) = i\/2FphysP^ 

TgiV^gt) = Uphys (46) 

The resulting lowest orders are 

Flo = F vlo = —BqF^ . (47) 
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Nc = 2 or pseudo-real representation Here we again need to be carefnl how we define 
the physical decay constant. We can choose to do using generators normalized to one using 
the original Dirac Fermions or generators normalized to one using the qi elements. 

With a Dirac fermion definition, the first Dirac (valence) quark labeled by 1 and the 
second by 2 , the definitions are 

{0|gi7/.75g2|M(p)) = i\/2FphysP;, 

{mi) = {imni + 1ri(1li) = Whys (48) 

M denotes a meson of that quark content with momentum p. The resulting lowest orders 
are 

Flo = vlo = —2BqF^ . (49) 

In terms of the g* the definitions are 

(0|gi7M?2 + Ql+Npl^^Q2+Np\M{p)) = 

'2{^i+Np,a^a,bCqiij Ql+Np ,a^abC Ql^h + Ql^a^abC qiJ^Np ,b 

7 Numerical examples and checks 

The main aim of this work is to provide the lattice work with the formulas and programs 
needed to do the extrapolation to zero mass. We therefore only present some representative 
numerical results. The numerical programs are included in the latest version of CHIRON, 

[Ml Eg. 

For the numbers presented we always use xi = 0.14^ GeV^, if not varied explicitly, and 
F = 0.0877 GeV as well as a subtraction scale p = 0.77 GeV. The length L for the finite 
volume has been chosen such that L x 0.14 GeV=3 or L 4.2 fm. 

The LEGs at NLO we choose to be those of the recent determination of [IH] with the 
extra LEG Lg = 0. The NNLO constants we have always put to zero. 

A number of numerical checks for the QGD case have been done. The unquenched 
infinite volume results for three flavours agree with the three flavour results of mm- 
The partially quenched results for masses and decay constants at infinite volume agree 
with the case dgea = l,dvai = 1 of [Ml [2S1 EZ]- The unquenched results for masses and 
decay constants at finite volume agree with [M] • The partially quenched results for masses 
and decay constants at finite volume agree with the case dsea = 1, dvai = 1 of [36] and 
finially the unquenched finite volume results for the vacuum expectation value agree with 
the results of [32] . 

In Fig. H] we show the mass squared for the infinite volume for all cases we have 
considered for three values of iV^. In general, as was already noticed in [TB] the cor¬ 
rections are larger for the larger values of Np. The corrections are also larger for the 
SU{2Nf) —>■ Sp{2Nf) case since this correspond to a twice as large number of fermions as 
the other cases. The partially quenched results shown in the right column are at a fixed 
value of Xi- That explains why the corrections do not vanish for 74 = 0. 


* V^FphysP fj. 
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The same types of results are shown for the decay constant in Fig. [21 The corrections 
are somewhat larger than for the masses but the convergence is typically somewhat better. 
The corrections for the vacuum expectation value shown in Fig. [3] are typically larger but 
with again a reasonable convergence from NLO to NNLO. 

We can now make similar plots for the hnite volume corrections. The overall size of 
them is as expected. The smallest mL is about two for the left hand sides of all plots. In 
the unquenched case the exponential falloff with the mass is clearly visible. The partially 
quenched cases contain a hxed mass scale Xi which is why the correction is more constant 
there, the stays at the mL = 3 point for the plots. The dips are caused by the finite volume 
corrections going through zero. The corrections to the mass are shown in Fig. IH the decay 
constant in Fig. |5]and the vacuum expectation value in FigO 

8 Conclusions 

We have calculated in the effective field theory for the three possible symmetry breaking 
patterns the NNLO order hnite volume and partial quenching effects to NNLO in the ex¬ 
pansion. The results satisfy a large number of checks agreeing analytically and numerically 
with earlier work that our results reduce to for some cases. The analytical part of this work 
relied heavily on FORM |19]. 

The analytical results are of reasonable length but given the total number of results we 
have included them as FORM output in a supplementary hie. They can also be downloaded 
from [5U] . 

The numerical programs have been included in CHIRON [38] version 0.54 which can 
downloaded from [39]. We have presented results in a number of cases with typical QCD 
values of the parameters. The results are of the expected sizes from earlier work in three 
havour ChPT. We hope these results will be useful for lattice studies of these alternative 
symmetry breaking patterns. 
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Figure 4: The absolute value of the finite volume correction to the physical mass squared 
divided by the lowest order mass squared for the unquenched (left) as a function of xi and 
the partially quenched case (right) as a function of Xa with xi = 0.14^ GeV^. Shown are 
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SU{Nf) SU{Nf). Middle line: SU{Nf) ^ SO{Nf). Bottom line: SU{2Nf) ^ 
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Figure 5: The absolute value of the hnite volume correction to the decay constant divided 
by the lowest order value Fq = Flo for the unquenched (left) as a function of xi ^md 
the partially quenched case (right) as a function of Xa with xi = 0.14^ GeV^. Shown are 
the NLO (p"^) and NNLO (p"^ + p®) results for three values of Np. Top line: SU{Nf) x 
SU{Nf) SU{Nf). Middle line: SU{Nf) ^ SO{Nf). Bottom line: SU{2Nf) ^ 
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value divided by the lowest order value vq = vlo for the unquenched (left) as a function of 
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